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Abstract Sun bear is found in Burma. Thailand, Vietnam, Sumatra, Borneo, 
North-east of India, Malaysia, and Yunnan and Sichuan of China and has been listed in 
1990 IUCN Red List of Threatened Animals as a vulnerable species. In this paper, a total of 
1174 bases of mitochondrial DNA sequence (397 bases of cytochrome b gene, 346 bases of 12S 
rRNA gene, 98 bases of tRNA genes, and 333 bases of D-loop region} from each of four sun 
bears from different collectuons have been sequenced. Comparing with the black rhino (RNA 
genes, our results suggest that the loop regions of tRNA genes evolye much fast than the stem 
regions. In the cytochrome b gene, 12S rRNA gene and tRNA gene regions, no sequence vari- 
ation among the individuals was observed, which indicates that genetic variability in the sun 
bear is relatively low. [mn the D-loop region. there are 13 and I sites. respectively, character- 
ized by transition and transversion among the 4 sun bears sequenced. The relationships among 
all the individuals were fully resolved based on mitochondrial D—loop sequences by using parsi- 
mony analysis. Our results suggest that the mitochondrial D-lcop region is very useful for popu- 


lation genetic study of the sun bear. 


Key words Sun bear, Mitochondrial DNA sequence. Conservation, Evolution 


There are seven species in the family Ursidae (Carnivora), among which sun bear 
{Helarctos malayanus) is the smallest, and inhabits the hot tropical forest. Sun bear is 
found in Burma, Thailand, Vietnam. Sumatra, Borneo. North-east of India, 
Malaysia. and Yunnan and Sichuan of China (Honacki e? al, 1982). It has been listed 
in 1990 IUCN Red List of Threatened Animals as a vulnerable species (IUCN. 1990) 
due to extensive destruction of the habitat and rapid decline of the population. 

Much effort should be made to save this endangered species. Population genetic ana- 


lyses is essential for a sound conservation management plan, and these questions should 
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be addressed; 1) How much genetic variability is there within populations? Is there sig- 
nificant difference among populations? 2} Are there gene flows among populations oc- 
curred recently? Do human activities have significant impact on the recent change of 
gene flows? 3) How big are the minimum viable population size (MVP) and the effective 
population size (Ne)? 

DNA sequences can provide abundant information on genetic structure of popula- 
tions (Nei, 1987; Avise, 1989). Rapid evolved mitochondria! DNA is a valuable mole- 
cule for investigating evolutionary relationships and genetic variations within and be- 
tween species (Wilson et al, 1985; Avise et al, 1987; Zhang et al, 1992, 1993a, b). 
Zhang ef al. (1991) found considerable divergence between mitochondrial! DNA haplo- 
types from two forms of Chinese pangolins, i. e. the brown and dusky. This suggests 
some genetic isolation between the brown and dusky Chinese pangolins. Studies on 
mitochondrial DNA in natural populations have shown that there is often concordance 
of phylogenetic relationship and geographic distribution of mitochondrial DNA 
haplotypes as a results of radiation or geographic barriers to gene flow (Lansman er al., 
1983; Zhang et al, 1993c). However, the phylogenetic relationships of mitochondrial 
DNA haplotypes in some species, such as house mice, do not always correspond to 
geographic distribution, which may be the result of human—mediated gene flow among 
the commensal! rodents (Ferris ef al., 1983). 

Wayne è: al. (1990) found highly divergent mitochondrial DNA haplotypes 
(p= 0.08} within a population of black—backals (Canis mesomelas}. They suggested that 
the great dispersal capability of carnivores might result in considerable gene flow and 
co-occurrence of divergent haplotypes within populations. 

Restriction enzyme analysis of mitochondrial DNAs from North American black 
bears, brown bears and polar bears suggests that there is the highest divergence within 
the American black bears and the lowest within the polar bears (Cronin et al. 1991). 
Mitochondrial DNA control region sequence analysis on the European brown bears re- 
veals a clear-cut difference between eastern European brown bears and the isolated popu- 
lations in the western Europe: more than 7% of the DNA sequence is different (Ta— ‘ 
berlet er al., 1992; Dorozynski, 1994). 

In our previous study, we analyzed mitochondrial DNA sequences of D—loop re- 
gion, cytochrome b, 12S rRNA, tRNA™™, and tRNA™ genes from all seven bear 
species and the giant panda. Our results suggest that the sun bear is the closest relative of 
the American black bear (Zhang et ail., 1993, 1994). 

The objective of this work was to investigate genetic variability within the sun 
bear, and to identify DNA region which is good for population genetic study of the sun 
bear. We present here the sequences from several mitochondrial DNA regions. 
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1 Materials and Methods 


1.1 Samples 

Four sun bear samples were collected. One was from Yunnan, China. Unfor-— 
tunately, the original localities of other three samples could not be determined. 
However, these three samples were from different collections. 
1.2 DNA Isolation and Nucleotide Sequencing 

Total genomic DNA was isolated from cultured fibroblasts, muscle, and liver tis- 
sues using the methods adapted from Lansman er a/. (1981). Conserved primer pairs 
L14724 + H15149 Urwin et al, 1991), L1091 + H1478, and L15926 + H00651 (Kocher 
et al, 1989) were used to amplify segments of cytochrome b, 12S rRNA genes and 
D-loop region (contained segments of tRNA?’and tRNA ™ genes) by PCR. PCR amplifi- 
cations were performed in a 100 wl reaction volume for 40 cycles using 2.5 units of Taq 
DNA polymerase (Perkin Elmer—Cetus) according to the manufacturer’s instructions, 
with denaturation for 40s at 95C, annealing for 1 min at 40-50CT, and extension for 
1-3 min at 72C . PCR reaction products were purified in 1.5% -2.0% low melting 
SeaPlaque agarose (FMC). The amplification product to be sequenced was cut out of the 
gel and melted in STE (0.1 M NaCl, 10 mM Tris-HCl, 1 mM EDTA, pH 8.0) at 
70€ for 15min, then extracted once with phenol and once with chloroform. DNA was 
precipitated with 2.5 volumes ethanol at ~20C for at least 2 hours and resuspended in 
water. Double stranded DNA was directly sequenced by the dideoxy method with 
Sequenase 2.0 (United States Biochemical) using heat denaturation as described by 
Zhang et al. (1993). Both the original PCR primers and internal primers could be used as 
sequencing primers. Products of sequencing reactions were electrophoresed on a 6% po- 
lyacrylamide / 50% urea gel. dried and subjected to autoradiography. The sequence was 
read using an IB] PC / GENE gel reader and sequence-—analysis software. All sequences 
were aligned using PC / GENE program version 6.6 (Intelligenetics, Inc.). 
1.3 Phylogenetic Analysis 

The variable D—loop nucleotide sequences were analyzed using the Phylogenetic Ana- 
lysis Using Parsimony (PAUP 3.1: Swofford. 1993). All nucleotide—substitution charac- 
ters were specified as unordered, and giving equal weight. American black bear was used 
for rooting the phylogenetic trees (Zhang ei al. 1994). 

Branch—and—bound searches were performed to ensure that all minimum~length 
trees were identified (Swofford, 1993). 

Bootstrap analyses consisted of 100 heuristic replications. As many as 100 trees were 
held for each bootstrap replication (Swofford, 1993). 


2 Results and Discussion 


A total of 1174 bases of mitochondrial DNA sequence (397 bases of cytochrome b 
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gene, 346 bases of 12S rRNA gene. 98 bases of tRNA genes, and 333 bases of D-loop 


region) from each sample have been sequenced. 


cytochrome b gene 


12S rRNA gene 
10 


Figure 上 


10 20 30 40 50 60 
ATGACCAACATC CGAAAAACC CACCCATTAGCTAAAATCATTAACAACTCACTTATTGAC 
MetThrAsniIleArgLysThrHisProLeuAlaLysI leIleasnAsnSerLeul leAsp 


70 80 90 100 110 120 
CTCCCAGC ACCATCAAACATCTCGGCGTGATGAAACTTCGGATCCC TCCTCGGAGTATGC 
LeuProAlaProSerAsnI leSerAlaTrpTrpAsnP heGlySerLeuLeuGlyvalcys 


130 140 150 160 170 180 
TTAATCCTACAGATTATGACAGGCCTATTTCTAGCCATACACTATACATCAGACACAACC 
LeurIleLeuGlnIleMetThrGlyLeuPheLeuAlaMetH isTyrThrSerAspThrThr 


190 200 210 220 230 240 
ACAGCCTTTTICATCAATCACTCATATCTGCCGAGACGTTCACTAC GGATGAATTATCCGA 
ThrAlaPheSerSerIleThrHisllecysArgas p¥ValHisTyrGlyTrplleIlearg 


250 260 270 280 290 300 
TATATACATGCAAACGGAGCC TCCATGTTCTTTATCTGCCTATTCATGCACGTAGGACGG 
TyxrMetHisAlaAsnGlyAlaSerMet PhePheI lecysLeuPheMetHisValGlyArg 


310 320 330 340 350 360 
GGTCTGTACTATGGCTCGTACCTATTCTCAGAAACATGRAACATCGGTATTATCCTCCTA 
GlyLeuTyrTy rGlySerTyrLeuPheSerGluThrTrpasnilec lyIleIleLeuLeu 


370 300 390 
TTTACAGTTATAGCCACCGCATTTATAGGATACGTCC 
PheThrValMetAlaThrAlaPheMetGlyTyrval 


20 30 40 50 60 
AAATAATTTATTAAACAARATTATTCGCCAGAGAACTACTAGCAACAGCTTAAAACTCAA 


70 BO 90 100 110 120 
AGGACTTGGC GGTGCTTTAGACCCCTCTAGAGGAGCCTGTICTGTAATCGATAAACCCCG 


130 140 150 160 170 184 
ATAGAC CTCACCACCTCTTGCTAATCCAGTCTATATACCGCCATCTTCAGCAAACCCTTA 


190 200 210 220 230 240 
AAAGGAACAAGAGTAAGCACAATCATCC CACA TAAAAAAGTTAGGTCAAGGTGTAATCCA 


250 260 270 280 290 300 
TGGGGTGGGRAGAAATGGGCTACATTTTCTTTTCAAGAACACTTACGAAAGTTTTTATGA 


310 320 330 340 
BRACTAAABACTAAAGGTGGATTTAGCAG TAAACCAAGAATAGAGAG 


Sequences of the sun bear mitachondmal! cytochrome b gene and [2S rRNA gene 


The predicted amino acid sequence of the cytochrome b protein is shown below in the three-letter anima 


acid code. There 1s no varauon among the four individuals studied in these two gene regions 
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19 20 30 49 50 
Sun bear I GTACATACTACTATTCCACCCAATGCCCTATTTATTTCATATATATCATT 
sun bear LI etme veces Neher etereeveasasees F Wala denier ew ew: We eleie eee c 
sun bear III a a ae ern C 
sun bear IV a 
American black bear ........C.T....9TT....TGOC.T.0.-5.Css L Ca 

60 70 B0 90 100 
sun bear I CTATGTACTGTACCATCGTAGTATGTITTCGAGTAC-TCCTTCTTTCATT 
sun bear II hh ee a moos 
Sun Dear IIT hh ee ee ee ol 
Sun bear Iv id ee ro me a a 0 


American black bear COR wate tere ree t sere rena eves Tole oo Tee Cones een. 


110 120 130 140 150 
sun bear 1 TITTTTCCCCCCTATGTACGTCETSCAPTAATGGCGTGCCCCATGCATAT 
fun bear II Worererer ror ooo oa a oe, 
sun bear IIL SAS es 
Sun bear IV CT 


American black bear te COs ccc ee rer eu ens ee te cnceseruccens 


169 179 1890 199 200 
sun bear I AAGCATGTACATATTATGCTTGGTCTTACATGAGGACCTACATTCCAAAA 
sun bear II 人 人 
sun bear III GT Ts 
sun bear TV Ts 


American black bear Ce 


210 220 239 240 256 
sun bear I GCTIGTITTAGGTGTATAGACTGTAAGCATGTATTTCACTTAGTCCGGGA 
Sun bear II Sr 
sun bear III "ev 
sun bear IV vevvevoeC。R.C.。.。。 a EN EEEE 00 6 wowed ee witin eae oe 8 
American black bear SR a 

260 270 280 2990 390 
sun bear I GCTIGATCACCAAGCCTCGAGAAACCAGCAACCCTTIGCGAGTACGTGTAC 
sun bear II Ca o me oie er ooo oe eeee 
Sun bear III er 
Sun bear IV CR i NELA AREE E E 
American black bear cs oo ss P 

a 310 320 330 

” sun bear I CTCTTCTCGCTCCGRGGCCCATGAAATGTGGGGG 
sun bear II 本 
sun bear III EEE ee ee ee E ee 
sun bear IV ne 


American black bear cso eres so oC 


Figure 2 Aligned mitochondria! D-loop sequences of the sun bear and American black bear 
Dots indicate identity with the sequence in the first row, and dashes indicate gaps inserted to 


maximize homologies. z 

2.1 Sequence variation in different mitochondrial DNA regions 
In the cytochrome b gene, 12S rRNA gene and tRNA gene regions, no sequence 
variation among the sun bear individuals was observed. Our results indicate that those 
gene regions can not be used as markers for population genetic study of the sun bear. 
The nucleotide sequences of mitochondrial cytochrome b gene and 12S rRNA gene, and 
the predicted amino acid sequence of the cytochrome b protein were presented in Figure 
1. The nucleotide sequences of mitochondrial (RNA? and tRNAT genes, and the pro- 
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posed cloverleaf structure of tRNA? were shown in Figure 3. 


tRNAThr ¢RNAPro 


10 20 30 40 50 
5'- AARARACGGAGAATATTCACCCTCCCCAAGACTCAAGGAAGAAGCAACAG 
60 70 ao 30 
CCCCGCTATCAGCACCCAAAGCTGATGTTCTATTTAAACTATTCCCTG -3' 
(a) 
tRNAPre 
A 
cG 
AU 
GU 
co 
co 
AU 
RU 
Pei CUUCG 
T ; 
A AyuuG IR NE: 
A II) nea 
U ,GAAC c 
A A 
ua U 
CG 
BU 
G-c 
CG 
U U 
U G 
UGG 
(b} 


Figure 3 Sequences of the sun bear mitochondrial tRNA genes 
There is no variation among the four individuals studied. “ + "indicates the position of the end and begining of 
tRNA7™ gene and tRNA* gene, respectively. A Watson-Crick base pair in a stem region is indicated by a dash. 
a} The partial sequence of tRNA?" gene and the complete sequence of t(RNA®™ gene, b) Proposed cloverleaf struc- 
ture of tRNA?” as inferred from gene sequence. 


In the D-loop region, there are 13 and | sites, respectively, characterized by tran- 
sition and transversion among the 4 sun bears sequenced. The sun bear D-loop se- ' 
quences have been aligned with that of the American black bear with introduction of one 
gap (Figure 2) (Zhang et al, 1994). 

Comparing cow mitochondrial tRNA genes with those of human, Anderson et 
ai.(1982) found that transition / transversion bias was especially marked in the stem re- 
gions of tRNA genes. They suggested that the constraints of certain structural require- 
ments favored transitions, since these allowed either Watson-Crick base—pairs to 
mutate to the other via a G-U or A-C intermediate. Since functional RNA genes, from 
the small tRNAs to much larger 18S and 26S rRNAs exhibit analogous secondary struc- 
tures that rely on base pairing between nucleotide positions. The stem regions of these 
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genes may be more conservative than the loop regions. For the sun bear tRNA? gene, 
there are 42 (64%) and 24 (36%) bases in the stem and loop regions, respectively. Com- 
paring with the black rhino tRNA*” gene (Jama er al. 1993), however. 70% of the 
substitutions (7 substitutions) occurred in the loop regions. Three transversions were ob- 
served, all of which occurred in the loop regions. Our results suggest that the loop re- 
gions of tRNA genes evolve much fast than the stem regions in the sun bear. 

2.2 Estimating phylogeny 

We choose the American black bear as the outgroup for rooting the phylogenetic 
tree. Even though the systematics of the ursine bear species has been the subject of con- 
siderable disagreement, recent sequence study suggests that the American black bear is 
the closest relatives of the sun bear (Zhang et al, 1994). Therefore, we consider that 
the American black bear is an ideal outgroup for resolving the relationships among the 
sun bear populations. 

Even through the original localities of three out of four samples could not be deter- 
mined, these samples were not from Yunnan and were from different sources. It is high- 
ly unlikely that they were from the same population. Therefore, our sample may repre- 
sent different sun bear populations. 

Phylogenetic tree was constructed based on the D-loop sequences with the parsimo- 
ny method using PAUP 3.1 (Swofford, 1993). A single shortest tree with 38 steps was 
identified (Figure 4). The consistency index and retention index are 0.95 and 0.87, 
respectively. All the branches were supported by > 70% bootstrap replications. The rela- 
tionships among all the individuals were fully resolved. Our results suggest that the 
mitochondrial D—loop region is very useful for population genetic study of the sun bear. 
2.3 Genetic variability 

It is not clear whether there is one or two subspecies in the sun bear, even though 
there is some body size difference between South~eastern Asia~Sumatra population and 
Borneo population. There is not any information on genetic variation between popula- 
tions. 

Mitochondrial cytochrome b gene sequence analysis has revealed considerable varia- 
tions among the populations in the Asiatic black bear, brown bear and polar bear 
(Zhang, unpublished: Zhang et al.. 1994; Zhang. 1995). In this study, however. no 
variation tn the mitochondrial cytochrome b gene and RNA gene regions between sun 
bear populations was observed. Our results indicate that genetic variability in the sun 
bear is relatively low. Interestingly, a few sequence variations among sun bear popula- 
tions were detected in the D-loop region, which suggests that different sun bear popula- 
tions belong to different maternal lineage. Possible explanation for the low genetic varia- 


bility could be a bottle~neck effect on sun bear populations in recently evolutionary histo- 


ry. 


OO 00D http://www.cqvip.com 


dng 动 物 学 研究 17 着 





2.4 Conservation relevance 


sun bear ] It is clear that saving a threatened 
3— sun bear I] 








species is much more complex than sim- 


sun bear H] : T k es 
ply setting aside habitat and raising a 
sun bear IY , 
23 en el bone few animals in the zoos (Woodruff, 
1989). Any sound conservation manage- 
Figure 4 Phylogenetic tree constructed based on ment plan should include the goal of pre- 


the D-loop sequences with the parsimony method 
using PAUP 3.1 (Swofford. 1993) 

A single shortest tree with 38 steps was identified. The 
consistency tndex and retention index are 0.95 and with low diversity are vulnerable to 


0.87. The lengths and bootstrap confidence values of changes in the environment. including n 
each branch are mdicated above and below the branch, 


serving as much as possible of genetic di- 


versity of the species, since populations 


resbechřel new pathogens and climatic changes re- 
sulting in diminished evolutionary poten- 
tial. Survival of the gene pool of endan- 
gered species will depend on a thorough understanding of the genetic variation of species 
and upon management strategies undertaken to protect and expand the surviving popula- 
tions. 

In our study, the genetic variability among sun bear populations studied is relative- 
ly low. Based on the mitochondrial DNA sequence as sole criterion. the populations in- 
vestigated could be managed as single conservation unit. Unfortunately, we could not de- 
termine if our samples represent the Borneo population. Further study with samples from 
the Borneo population and other geographic populations will provide useful information 
on the conservation of the sun bear. 
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BE 马 来 能 在 IUCN 红皮书 中 被 列 为 受 胁 动物 ， 其 保护 受到 广泛 的 关注 。 本 文 研究 4 
只 马 来 能 的 线粒体 DNA 序列 ， 其 中 1 只 来 自 云 南 ， 其 余 3 只 产地 不 详 ， 但 来 自 不 同 的 搜集 渠 
道 。 对 于 每 个 个 体 ， 我 们 测定 了 397 bp 的 细胞 色素 b 基因 、346bp 的 12S rRNA 基因 、 
98 bp BY tRNA 基因 和 333 bp 的 D 环 区 序列 ， 共 计 1174bp。 经 与 黑 犀 牛 序列 比较 ， 发 现 
RNA 基因 的 空间 结构 对 基因 的 进化 有 显著 影响 ， 环 区 的 进化 明显 快 于 蔡 区 的 。 对 于 细胞 色素 
b 基因 、12S rRNA 基因 和 RNA 基因 ， 在 马 来 熊 个 栖 间 未 发 现 序 列 变异 。 这 一 结果 提示 ， 马 
来 艇 群体 的 遗传 变异 程度 较 恢 。 在 D 环 区 ， 有 13 和 1 个 位 点 分 别 出 现 转 换 和 其 换 。 根 据 D 
环 区 序列 ， 我 们 采用 简约 法 确定 了 马 来 能 群体 问 的 进化 关系 。 我 们 的 结果 表明 ， 钱 粒 体 DNA 
的 DD 环 区 是 研究 马 来 能 群体 遗传 结构 十 分 有 效 的 遗传 标记 。 
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